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The first lattice QCD result on the nuclear force (the NN potential) is presented in the quenched 
level. The standard Wilson gauge action and the standard Wilson quark action are employed 
on the lattice of the size 16 3 x 24 with the gauge coupling j3 = 5.7 and the hopping parameter 
K = 0.1665. To obtain the NN potential, we adopt a method recently proposed by CP-PACS 
collaboration to study the %% scattering phase shift. It turns out that this method provides the 
NN potentials which are faithful to those obtained in the analysis of NN scattering data. By 
identifying the equal-time Bethe-Salpeter wave function with the Schrodinger wave function for 
the two nucleon system, the NN potential is reconstructed so that the wave function satisfies the 
time-independent Schrodinger equation. In this report, we restrict ourselves to the J p = + and 
1=1 channel, which enables us to pick up unambiguously the "central" NN potential V C entrai('")- 
The resulting potential is seen to posses a clear repulsive core of about 500 MeV at short distance 
(r < 0.5 fm). Although the attraction in the intermediate and long distance regions is still missing 
in the present lattice set-up, our method is appeared to be quite promising in reconstructing the 
NN potential with lattice QCD. 
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1. Introduction 

The nuclear force (the NN potential) serves as one of the most important building blocks in 
nuclear physics. Its attractive part in the intermediate distance region plays an essential role to 
bind nucleons in nuclei. On the other hand, its strong repulsive core at short distance ensures the 
stability of heavy nuclei by preventing them from collapsing and leads to the celebrated saturation 
phenomena of nuclear matter. Importance of the nuclear force goes even beyond the nuclear struc- 
ture: It can have a great influence on the physics of compact stars and the supernova explosion 
through the equation of state at high baryon densities. 

So far, enormous effort has been devoted to understand the nuclear force theoretically [jj]]. 
The long distance region (r > 2 fm) of the nuclear force is dominated by the exchange of a single 
pion, i.e., the lightest elementary excitation in the QCD spectrum. The intermediate distance region 
(1 ~ r !& 2 fm) receives a significant contribution from multi-pion exchanges and from heavy meson 
exchanges such as p, ft), and "a". The understanding of the short distance region (r ;$ 1 fm) is most 
retarded, but it is believed to be intimately related to the quark-gluon structure of the nucleons. 

Hence, the first principle QCD calculation of the nuclear force (in particular the short distance 
region) has been desired for a long time. Recently, a method attempting to obtain the origin of 
the repulsive core on the lattice was proposed, where the nucleon is assumed to be composed of 
heavy-light-light quarks so that one can define the relative distance between the nucleons in terms 
of the distance between the heavy quarks [g]. (The same technique was also used in Ref. [||].) In 
this report, we introduce a totally different approach to the NN potential on the basis of a method 
recently proposed by CP-PACS collaboration in the context of Tin scattering phase shift For 
simplicity, we confine ourselves to the channel J p = + and 7 = 1, which makes us possible to 
obtain unambiguously the so-called central potential Ventral (*")• First, we construct the equal-time 
Bethe-Salpeter (BS) wave function for the NN system. It is then identified with the non-relativistic 
Schrodinger wave function for two nucleons interacting at low energies. V cen trai( r ) is reconstructed 
so that the wave function satisfies the non-relativistic Schrodinger equation. In this way, we do not 
have to introduce heavy quarks to define the relative distance unlike the method proposed before. 

The reconstructed Ventral M has a repulsive core of about 0.5 GeV in the short distance region 
(r < 0.5 fm). On the other hand, the attraction in the intermediate and long distance regions is still 
missing, which may be attributed to combined artifacts of the heavy pion mass (m n ~ 0.5 GeV), 
the small spatial volume (L ~ 2.2 fm) and statistics. 

2. The formalism 

We consider the equal-time Bethe-Salpeter (BS) wave function <p(x;k) for NN system in the 
J p = + and 1=1 channel as 

p a {x) = ea bc {u T a Cy 5 d b )u Ci a, 
np(x) = £abc{u T a Cy 5 db)d c ^, 

where a, b and c denote the color indices, a and /3 the Dirac indices, and C = 74/2 the charge 
conjugation matrix. p a {x) and np(y) denote composite bispinor fields for proton and neutron, 



p a (R-x+X)n p (X) 



pn;k 



(2.1) 
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which are represented with Dirac's convention, x plays the role of the spatial separation of the 
proton from the neutron. The summation over R is performed for the cubic transformation group 
O to pick up the s-wave orbital contribution of NN system. Note that, due to the parity selection 
rule, only s-wave contribution is allowed for J p = + . The spinors of the proton and the neutron 
are combined into scalar in the non-relativistic manner with {02) a p- Since we are interested in 
the significantly low energy region, the lower components of the nucleon bispinors, i.e., p a and 
nn, are negligible. Since this wave function is symmetric in orbit and anti-symmetric in spin 
configuration, it has to be symmetric in iso-spin space due to the Pauli statistics, i.e., this wave 
function is iso- vector. The summation over X is performed to project the total spatial momentum 
to zero, i.e., P = 0. The Liischer's k is introduced as a relative momentum outside the range of the 
NN interaction, playing the role of "asymptotic relative momentum" between the proton and the 
neutron in the scattering theory in the infinite volume limit. 

The BS wave function Eq. ( |Q| ) is obtained, in the lattice QCD formulation, from the large t 
behavior of the four point nucleon correlator as 

'o\ (2.2) 



F NN (x,y,t;to) = (0 p a (x,t)np(y,t)p a >(0,to)np>(p,t ) 
= ^(0\p a (x)n p (y)\n) A n e - E ^\ 

n 

where E n denotes the energy of the state \n), to the time-slice on which the source is located. 



A„ = (n\p a {0)np(0)\0) plays the role of an "overlap". Eq. (2T) satisfies the Bethe-Salpeter (BS) 
equation, since it is obtained as a solution to the BS equation for the NN system in the BS frame- 
work. By using the procedure of the non-relativistic reduction, the BS equation reduces to the 
effective Schrodinger equation |Q] as 

(V 2 + k 2 )<p(x;k) =m N J d 3 yV k (x-y)<p(f,k), (2.3) 

where the nucleon mass m^ is introduced in the r.h.s. for later convenience. The interaction kernel 
Vic(x — y), in general, depends on k, which is the main reason why Eq. ( p73| ) is referred to as "effec- 
tive" Schrodinger equation. At low energies, Vk can be approximated as Vk{x — y) ~ Vk=o{x — y). 
We are left with the Schrodinger equation for the NN system as 

-—V 2 $(x;k) + J d 3 yV NN (x-y)(t>(f,k)=E$(x;k), (2.4) 

where /i = m^/2 denotes the reduced mass of the nucleon, and E = k 2 /(2ju) denotes the non- 
relativistic "energy". 

According to the standard nuclear physics [S], the NN interaction Vnn is parameterized in the 
low energy region as 

V NN = Vcentral(r) + Kmsor W § 12 + V LS (r)L ■ S+ 0(p 2 ), (2.5) 
S\2 = 3 {O p - O n ), 

L = -ixxV, S={dp + 6 n )/2, 
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where r = \x\ denotes the distance between the proton and the neutron. a p and a„ denote the spin 
Pauli matrices for the proton and the neutron, respectively. L represents the relative angular mo- 
mentum operator of the proton from the neutron. Ventral (r) is referred to as the "central" NN force, 
which only depends on the distance r. It is considered to be the most important local interaction in 
the NN interaction. (Interactions which do not involve any derivatives are referred to as "local" .) 
VtensorM is referred to as the "tensor" force, which is the second important local interaction. VlsW 
is referred to as the "spin orbit interaction" , which is the most important non-local term. Finally, 
0{p 2 ) represents remaining non-local contributions, which is often considered to be less important. 

Restriction to the J p = + channel permits a further simplification. Because only the s-wave 
orbital component is allowed in this quantum number due to the parity selection rule, the contribu- 
tions from the second and the third terms vanish. Since, then, V^n receives the contribution only 
from the central term as 

VnN - VcentralO"), (2.6) 



Eq. ( |2.4| ) reduces to 

V 2 

-—<l>(x;k)+V ceDtX ti(r)<l>(x;k)=E<i>(x-,k). (2.7) 
Since Kentrai (?) is a simple multiplication operator, we can rearrange it in the following way: 

1 W 2 <b(x\k) 

Y ~M =E+ 2t-m>- (2 ' 8) 

This relation states that, by identifying the BS wave function for NN, which is obtained from 
lattice QCD calculations, as the Schrodinger wave function, Ventral (r) can be reconstructed so that 
the Schrodinger equation is satisfied. It should be emphasized that, by construction, nothing is 
required on the quark mass to define the fixed separation between the two nucleons. Since the NN 
potential is reconstructed from the NN wave function, this method is expected to provide such NN 
potentials that are faithful to those obtained in the analysis of the NN scattering data. 



3. The numerical calculation 

We employ the standard Wilson gauge action at the gauge coupling j3 = 5.7 on the 16 3 x 24 
lattice together with the standard Wilson quark action with the hopping parameter fc = 0.1665. 
The lattice spacing is determined from p meson mass in the chiral limit as a -1 = 1.44(2) GeV 
(a ~ 0.137 fm). The physical volume of our lattice is (2.2fm) 3 x (3.3fm). These parameters 
reproduces m K ~ 0.53 GeV, m p ~ 0.88 GeV, and — 1 .3 GeV. (See Ref. [0] for detail.) The global 
heat-bath algorithm is used to generate the gauge configurations. After skipping 5000 sweeps 
for thermalization, the gauge configurations are picked up every 100 sweeps. Totally, 160 gauge 
configurations are used for the measurement. Dirichlet boundary condition is imposed on the quark 
fields along the temporal direction on the time-slice t = 0. To enhance the ground state contribution 
of the NN wave function, we adopt a spatially extended source with the gaussian smearing method 
with the size p ~ 0.4 fm on the time-slice t = 5. The BS wave function for NN is measured on 
the time-slice t = 10, which is determined from a plateau appearing in the NN effective mass plot 
shown in Fig. (|l|). We keep in mind that the ground state saturation of the BS wave function is quite 
important. 
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Figure 1: The effective mass plot. Dots denote the effective mass of the NN system (J p 
whereas crosses denote twice the effective mass of the nucleon. 



0+, /= 1), 



2 . 



o 

■H 

+J 
O 

c 

m 

<u 
> 

3 



1 . 5 



1 . 



. 5 




mi 1 1 



. 



. 



. 5 



1 . 

r [fm] 



1 . 5 



2 . 



Figure 2: The lattice QCD result of the NN wave function (J p = 0+, / = 1). 
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Figure 3: The lattice QCD result of the NN potential V C entrai(>")- 



Fig. (g) shows the lattice QCD result of the BS wave function for NN. It is normalized with 
its value at the origin 7 = 0. Not only the on-axis data , but also the off-axis data are plotted 
all together, due to which the horizontal axis extends beyond L/2 ~ 1.1 fm to V3L/2 ~ 1.9 fm 
(L ~ 2.2 fm is the spatial lattice extension). We observe that the wave function shrinks near the 
origin suggesting the existence of repulsion at short distance. 

Fig. (||) shows the lattice QCD result of the NN potential V ce ntrai (r) for J p = 0+ , / = 1 channel. 
The simplest nearest neighbor numerical Laplacian for Eq. ( |2.8[ ) is employed in reconstructing 



VcenterM- The zero adjustment due to E = j^k 2 in Eq. (18) has not yet been performed because 
of the large noise. In principle, Liischer's k can be obtained from the difference of the effective 
masses, i.e., »iAw, e ff — ^niN,eS (cf. Fig. (|])), or the NN wave function in the large r region [ph. So 
far, the data is too noisy to obtain it. We observe that, corresponding to the shrink in Fig. ([J), 
^central ( r ) nas a clear repulsive core of about 500 MeV in the short distance region r & 0.5 fm. On 
the other hand, we do not find a significant attraction in the intermediate distance region, which 
may be attributed to (i) the still large noise, (ii) finite volume effect (L ~ 2.2 fm is too small for 
NN), (iii) the large pion mass (m n ~ 0.53 GeV). In order to gain physics insights into the origin of 
the repulsive core at short distance, we need more information such as the quark mass dependence, 
the flavor dependence, etc. 



4. Summary and discussion 

To study the nuclear force, we have applied a new method recently proposed by CP-PACS 
collaboration in the context of the %% scattering phase shift. In this method, constraint on the quark 
mass to define the distance between the two nucleons is not necessary. Because the NN potential 
is reconstructed from the NN wave function, it is expected to provide NN potentials, which are 
faithful to those obtained from the analysis of the NN scattering data. By restricting ourselves to 
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the quantum number J p = + and 7 = 1 , we have reconstructed the central part of the NN-potential 
V cen trai('") fr° m the Bethe-Salpeter wave function for NN obtained by using lattice QCD. We have 
seen that Ventral ( r ) has a clear repulsive core of about 500 MeV in the short distance region r ;$ 0.5 
fm. However, we have not found a significant attraction in the intermediate distance region, which 
may be attributed to the poor statistics, the small lattice volume (L ~ 2.2 fm) and the large pion mass 
m % ~ 0.53 GeV. We are currently performing a lattice QCD calculation for V cen trai ( r ) using a better 
statistics on a larger lattice volume with a lighter pion mass, which will be presented elsewhere. 
Our method have appeared to be quite promising in reproducing the NN potential with lattice QCD. 
It is interesting to apply our method to other channels to construct Vus(r) and Vtensor(/)- Also, it 
is important to apply our method to hyperon-nucleon and hyperon-hyperon forces, where only a 
limited number of experimental data are available so far. 
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